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We show that the perturbative expansion of the two- level correlation function, R(uj), in disordered 
conductors can be understood semiclassically in terms of self-intersecting particle trajectories. This 
requires the extension of the standard diagonal approximation to include pairs of paths which are 
non-identical but have almost identical action. The number of diagrams thus produced is much 
smaller than in a standard field-theoretical approach. We show that such a simplification occurs 
because R{u>) has a natural representation as the second derivative of free energy F{uj). We calculate 
R{u3) to 3-loop order, and verify a one-parameter scaling hypothesis for it in 2d. We discuss the 
possibility of applying our "weak diagonal approximation" to generic chaotic systems. 



PACS numbers; 73.23.-b, 73.20.Fz, 05.45-Hb 

The relationship between the quantum properties of 
disordered and classically chaotiy systems has been the 
focus of much recent researchETO. The main reason for 
this is that both types of system show the same un- 
derlying behavior - their energy spectra in appropriate 
regionS|-|liave statistics given by random matrix theory 
(RMT)lI'EI. The spectral properties in these universal 
regions are well understood. The challenge now is to 
understand whether generic features also emerge in the 
deviations from universality. Answering this question is 
made difficult because two very different languages are 
used to characterize the two types of system. In quan- 
tum disordered systems one averages over all realizations 
of disorder to obtain_an effective field theory: the non- 
linear sigma modelo^E3. In quantum chaotic systems, (i.e. 
quantum systems which are chaotic in the classical limit) , 
one considerS|-a particular system using the Gutzwiller 
trace formulaEil, which involves a sum over classical pe- 
riodic orbits. Each language has its own strengths and 
weaknesses. The field theory technique is rather formal 
and somewhat opaque physically, but has a well-defined 
perturbation expansion. The trace formula appears to 
be more physically transparent in that one is summing 
over classical trajectories. However, this sum is difficult 
to perform, its convergence properties are not well under- 
stood, and no controlled expansion is currently available. 
It is therefore natural that one should attempt to use the 
strengths of one language to compensate the weaknesses 
of the other. 

One recent attempt was to try to carry over the power- 
ful calculational techniques dejweloped for the disordered 
systems to the chaotic systemsoO. In the most consistent 
form, one averages over a certain energy interval to gener- 
ate a field-theoretical functional in the Wigner represen- 
tation (the "ballistic" sigma-model). From this one can 
derive, e.g., the two- level correlation function for quan- 
tum chaotic systems. In the universal ergodic regime, the 
non-peiturbative derivation is equivalent to one which 
EfetovH developed for disordered systems, as expected. 



In addition the leading-order perturbative termfl in the 
non-ergodic regime gives a result equivalent to that ob- 
tained in the diagonal approximatioiii3 to the Gutzwiller 
trace formulaO. The same, approximation when applied 
to the disordered systemad also gives the leading order 
perturbative term first derived in Ref. ^ One therefore 
expects that higher order perturbative terms in quantum 
chaotic systems would be analogous to those occurring 
in disordered systems (where they are known as weak- 
localization corrections). The ballistic sigma model is, 
however, ill-defined in the ultraviolet limit, making per- 
turbative analysis currently ambiguous. 

On the other hand, various authors&0§ have pro- 
ceeded in the opposite direction, and developed semiclas- 
sical methods for the quantum disordered system. The 
goal is to gain a more intuitive picture of phenomena such 
as weak localization and universal conductance fluctua- 
tions. 

In this paper we show how the diagrammatics for spec- 
tral correlations in disordered conductors can be rewrit- 
ten in terms of particle trajectories which self-intersect 
in real space. Each diagram for the two-level corre- 
lation function is represented as consisting of two tra- 
jectories that are identical everywhere except at self- 
intersection regions, where they are rejoined in different 
ways. As the two trajectories are identical for most of 
their length, they are phase coherent, and interfere con- 
structively. The perturbation parameter, 1/g, where g 
is the the dimensionless conductance, can be understood 
as the probability for a self-intersection to occur. The 
great advantage of this approach is the drastic reduction 
in the number of diagrams which occur in a given or- 
der of perturbation theory This reduction arises because 
our approach does not distinguish the starting points of 
the (closed) trajectories, whereas the standard approach 
does. We show that when translated into field theory lan- 
guage, this means that the two-level correlation function, 
i?(cij), is the second derivative with respect to energy, uj, 
of the free energy, F{uj). As a particular problem we 
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calculate the two-level correlation function to three-loop 
order in perturbation theory, which yields the leading or- 
der level statistics for the unitary system in 2d. Finally 
we show that the three-loop results for R{lu) in both uni- 
tary and orthogonal systems can be derived from a one- 
parameter scaling hypothesis in which the rcnormalizcd 
conductance g{i-j) is substituted into the one-loop result. 
Note that in our approach we are using the trajectory pic- 
ture to classify the field theory diagrammatics in contrast 
to the work discussed in the previous paragraph where 
field theory is being used to classify periodic orbits. 

We consider the relation of the picture of trajectories 
to the perturbation expansion of R{uj) which is defined 

by 

R{oj)^ ^(^u{E + Lj)i^{E)'j -1. (1) 

Here I'iE) is the density of states per unit volume which 
can be written, given the spectrum {En} of the system, 
as 

n 

and I' = {v{E)) = 1/{AL'^), where A is the mean level 
spacing. Here (. . .} denotes the averaging, either over all 
realizations of disorder for a disordered system or over a 
certain energy window for a chaotic system. Such energy 
windows are always narrow enough so that v has no en- 
ergy dependence. The semiclassical approach allows one 
to write v{E) as a sum overj^assical paths by means of 
the Gutzwiller trace formulate, 

iy{E) = ME) exp {^Sp{E)/h) (3) 
p 

where Ap{E) and Sp{E) are the amplitude and action of 
the p-th periodic orbit at energy E. One can then substi- 
tute this semiclassical formula into the definition of the 
two- level correlation function, R{uj), and take the Fourier 
transform to obtain the spectral form factor, K{t), 

K{t)= [ ^e^^{-iujtlh)R{u). (4) 
J 2nn 

Since all factors vary slowly with E except the action 
Sp{E), one expands Sp{E + oj) to first order to obtain 

K{t) = Y.ApA;^^''''~''^'>/H[t~^{Tp + T,)] (5) 

P,9 

where Tp{E) = dSp{E)/dELl£ the period of the orbit p. 

To proceed further BerryEj introduced the diagonal ap- 
proximation assuming that only terms with p — q con- 
tribute to the sum in Eq. (|). The reasoning is that 
terms with p ^ q have randomly varying phases and can- 
cel each other. In the disordered metal such a diagonal 
approximation was shownlJ to-4'eproduce the leading or- 
der perturbative contribution!!^ to R{uj). 



The question that we address here is how to get reg- 
ular corrections to the diagonal approximation which 
are due to contributions of relatively long trajectories. 
Such corrections are of fundamental importance, espe- 
cially for two-dimensional disordered systems: the lead- 
ing order term corresponding to the diagonal approxi- 
mation vanishescZI for d = 2. However, the possibility 
of getting the regular corrections could be interesting in 
a much wider context. To evaluate the corrections one 
has to consider pairs of paths with almost equal actions, 
Sp w Sq, like in the diagonal approximation, but allow 

p^q- 

To construct the weak diagonal approximation which 
allows for such pairs of trajectories we start with a path 
that self-intersects in real space at one or more points. 
At a point of intersection we break up the path and make 
different trajectories by joining the pieces together in dif- 
ferent ways. We see at once that this gives rise to a per- 
turbation expansion, with the perturbative order given 
by the number of loops in real space created by self- 
intersections, and the perturbation parameter given by 
the probability of having a self-intersection. In the case 
of a disordered electronic system, we will show that this 
expansion is none other than the usual field-theoretical 
loop expansion. Before we do this, let us speculate on 
the nature of this expansion for a generic chaotic system. 

In this case the above picture of self-intersecting trajec- 
tories requires clarification. First of all, the picture of tra- 
jectories is formulated in phase space rather than in real 
space. Of course, classical trajectories which are identical 
along part of their length in phase space cannot diverge 
in the way discussed above ~ this process is quantum me- 
chanical in nature. The uncertainty principle means that 
phase space is coarse-grained into boxes with size of order 
h'^. Two trajectories which were originally nearly iden- 
tical (i.e. passing through the same phase space boxes) 
start to deviate significantly at the Ehrenfest time, ^Ehr, 
due to quantum processes such as diffraction. Then we 
can make pairs of piecewise identical trajectories as de- 
scribed above. For this picture to make sense we need 
to work at time-scales greater than ^Ehr- This has been 
noted in Ref. |l^ which recently considered the weak lo- 
calization correction in a disordered system where the 
Ehrenfest time is determined by the diffraction at ran- 
domly distributed scattering spheres. In the usual model 
of a disordered conductor with point-like scatterers ^Ehr 
coincides with the elastic scattering time r since the par- 
ticle's direction is totally randomized after each scatter- 
ing. In general, we expect that the classification of tra- 
jectories in the weak diagonal approximation might be 
valid for chaotic system with t-Ehr ^ iorg where torg is 
the ergodic time scale at which the trajectory samples 
the entire phase space. This condition is necessary in 
order to treat the self-intersection region as a perturba- 
tion. It simply means that the length of almost identical 
regions of the two trajectories is much larger than the 
length of the self-intersection region. 
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At time scales t ^ tag, one expects completely uni- 
versal behavior which is described by the random matrix 
theory.lla This has been conjectured bv-jBohigas et alQ 
and partially proved by Andreev et alH'tl This univer- 
sal behavior is completely non-perturbative and cannot 
be described by the diagonal approximation. Moreover, 
even for t ^ tn, where tn = ?i/A is the Heisenberg time, 
the diagonal approximation does not reproduce the de- 
viation from the linear behavior of K{t) known from the 
random matrix theory. We will show that we can ex- 
actly allow for this deviation within our weak diagonal 
approximation. 

At time scales tshr ^ i ^ tcrgj one expects non- 
universal corrections to the diagonal approximation. It 
is not clear whether this is true for a generic chaotic sys- 
tem. Attempta23 to generalelhese corrections within the 
ballistic nonlinear a modelM~El have led to results vanish- 
ing to all orders in perturbation theory. These results 
are dubious, however, because they depend strongly on 
the short-time regularization procedure. In particular, 
there exist regularizations where the corrections do not 
vanish and are totally governed by short-time processes 
at time scales t ^ tEhr- In addition, the ballistic a model 
itself might not take into account all relevant quantum 
processes. 

In the disordered conductor these corrections can be 
calculated both for t ^ terg and for t <ti icrg- Here 
tcrg = / D is the time taken for an electron to diffuse to 
the boundaries. We will show that the trajectory picture 
is naturally related to a standard diagrammatic expan- 
sion and, moreover, allows considerable simplifications in 
this expansion. 

We will now draw all closed trajectories with up to 
two self-intersections regions and show their relation to 
the diagrams up to three-loop order. Let us first rewrite a 
pair of identical (up to discrete symmetries of the ensem- 
ble) trajectories-pwhich correspond to Berry's diagonal 
approximationJ13id as a diagram in field theory language. 




(a) (b) 

FIG. 0. (a) The 1-loop diagram for the free energy F(lli), 
which consists of pairs of trajectories that are identical up to 
time reversal symmetry. In the field theory language it con- 
sists of a single closed wavy line, (b) The same diagram in 
the disordered metal. Since the closed wavy line in (a) con- 
sists of a sum over impurity ladders, and the diagram with n 
impurity ladders has symmetry factor 1/n, the ladder gives 
logarithmic contribution — ln(Z)g^ — iuj). 




FIG. y. The 2-loop diagram for the free energy F{uj) cor- 
responding to paths with one self-intersection. First we show 
the two different ways of linking up the paths at the self- 
intersection point. Then we rewrite this in the field theory 
language so that regions where the two paths are identical 
become wavy lines, whilst the region of self-intersection be- 
comes a closed box. Finally we twist the latter diagram into 
its usual form. 

The pair of identical or time-reversed trajectories be- 
comes a diffusion or Cooperon propagator, respectively, 
drawn as a wavy line. This yields the one-loop diagram 
shown in Fig. la. 

The two-loop case is shown in Fig. 2 and involves one 
self-intersection region. At this point the paths can be 
linked up in two different ways, and the diagram con- 
sists of these two trajectories. Again we rewrite in field 
theory language with the identical or time-reversed por- 
tions of path becoming diffusion or Cooperon propaga- 
tors. The region around the self-intersection becomes an 
effective interaction between these propagators known as 
a Hikami boxEII. If we put arrows on the trajectories to 
show their direction of traversal we see that they go in 
the same direction over part of the diagram. This means 
that time-reversal invariance is required for their actions 
to be identical. Hence there is no two-loop term in the 
unitary case, where time-reversal invariance is broken; 
in fact, all even loop contributions vanish for the same 
reason. 

The three-loop case is shown in Fig. 3, and involves 
two self-intersections. The situation here is more compli- 
cated since there are three distinct topological ways for 
two self-intersections to occur. They can occur either at 
the same point, or at two different points A and B. In 
the latter case there are two distinct orderings in which 
the trajectory can be traversed - AABB and ABAB. 
Finally there is more than one way of linking up the par- 
tial paths to form trajectories. There are five diagrams 
in total: F^a has AABB form; F^i, and F^c have inter- 
section at only one point; F^d and F^^ have ABAB form. 
Putting arrows on the trajectories we find that only i^sb 
and F^d contribute to the unitary case. 

We see that the above procedure can obviously be gen- 
eralized to any order of perturbation theory, and is a 
powerful way of ensuring that all contributions have been 
considered. 

Let us now relate these pictures of semiclassicaLtra- 
jectories to the standard diagrammatic approachll3 for 
evaluating R{uj). The starting point is the expression for 
the density of states, Eq. (|2|), in terms of electron Green's 
functions. 
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FIG. |3| The five 3-loop diagrams for the free energy F(uj) corresponding to paths with two self-intersections. We follow the 
format of Fig. 2, showing first the pairs of distinct trajectories; then rewriting these with wavy lines representing regions where 
the two paths are identical, and closed boxes at the self-intersection points; then finally twisting the latter diagrams in their 
usual form. 
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J d-^r [G^(r, r; E) - G'^(r, r; E)] 



It follows from this expression and the definition (|l|) that 

R{io) = ^^eJd^rd''r'{G''ir,r;Ei)G^{r',r';E2))^ 



uj = El- E2 



(6) 



(7) 



Here only the connected average of G^G^ remains: 
the sum of all unconnected averages is absorbed by —1 
in the definition (|l|), while the connected average of 
G^G^ + G^G^ vanishes. 

Equation (7) directly corresponds to the semiclassical 
expressions (i) and since G^'{r,r;E) is the quan- 
tum mechanical propagator for a particle of energy E to 
start and finish at the same point r. We can then rep- 
resent G^{r,r;E) as the sum of all closed paths while 
G^{r, r; E) will be the sum of all closed paths traversed 
in the opposite direction. Averaging over disorder results 
in the vanishing of all contributions apart from those 
corresponding to pairs of coherent paths. This coher- 
ence arises when the pair of paths involves scattering 
oS the same impurities although not necessarily in the 
same order. Pictorially, scattering off the same impuri- 
ties is shown by impurity lines between G^ and G^. In 
the lowest perturbative order, the closed paths for G^ 
and G^ are identical, apart from their different starting 
points r and r', while the impurity lines form a ladder 
which corresponds to a diffusion propagator 
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< 1, WT < 1 (8) 



or its time- reversed counterpart, a Cooperon. 

The key feature of the diagram for generated in 

such a way. Fig. 4, is that the starting points of the two 





FIG. § The field theory diagram for the lowest order con- 
tribution to the two-level correlation function R{lo) in the 
standard approach. This is identical to Fig. 1 except that 
the starting points of the two trajectories, r and r' are dis- 
tinguished. In general inserting points in such a way yields 
many more diagrams for R{io) than for F(ll>). 



paths, r and r', are fixed in the beginning, even though 
we finally integrate over all r and r' in Eq. (^. Con- 
trarily, there are no starting points on the semiclassical 
diagrams of Fig. 1. Obviously, one can obtain the dia- 
gram of Fig. 4 from the semiclassical picture of Fig. 1 by 
inserting points r and r' into different electron lines in all 
possible ways. Inserting an external point is equivalent to 
taking a derivative with respect to energy: inserting r re- 
places G^{Ei) by G^{Ei)G^{Ei), which can be achieved 
by the action of —d/dEi. Similarly, inserting r' gives the 
same result as —djdEi. After averaging R(uj) depends 
only upon lo — E\ — E2. Therefore, one may represent it 
as the second derivative of a certain function F{uj), 



R{uj) 



(9) 



We will show below that this F{uj) corresponds to the 
'free energy' of the appropriate field-theoretical func- 
tional. Note that such a representation of R{uj) has previ- 
ously been used in thej-evaluation of first-order diagrams 
in the ballistic regimeE3. Here we show that Eq. (||) is 
valid to all orders of perturbation theory. 

Higher order perturbative contributions to arise 
from two paths where scattering occurs from the same 
impurities but in a different order. The parts of the di- 
agrams where the sequence of impurity scatterings coin- 
cides (or is time reversed) for the two paths are again 
represented by diffuson (or Cooperon) ladders. These 
ladders are connected by Hikami boxes which represent 
the change in the scattering sequence. In terms of paths, 
this corresponds to the regions of self-intersection de- 
scribed above. The only difference between these higher 
order contributions to Riuj) and semiclassical diagrams 
in Figs. 2 and 3 is the necessity to distinguish starting 
points r and r'. Therefore these semiclassical diagrams 
describe higher order contributions to F(lu), and R{uj) 
can be obtained with the help of Eq. (j^. 

Finally let us further justify the above pictorial discus- 
sion with a rigorojis drrivation using the standard field- 
theory machinerylj'HElrEj. We perform the average over 
disorder by the replica trick. This can be done with ei- 
ther bosonic (commuting) or fermionic (anticommuting) 
variables, and both must yield the same results for physi- 
cal quantities. In this paper we will use bosonic variables. 
The two-level correlation function is then given by 



Ri^) 



lim — 



VQ / d'^rTr[AQ(r)] exp {-T[Q;uj]) - 1 



(10) 
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where is the non-Hnear sigma model functional 

T[Q-uj] J '^''''^^ [^(^Q)' - (11) 

The matrix Q has rank 2N: we made N replicas each for 
retarded and advanced Green's functions. It satisfies the 
standard saddle-point conditions: 



I . 



Tr(Q) = . 



The matrix A is diagonal with elements +1 for retarded 
indices and —1 for advanced indices. Then it follows that 

Tr(AQ(r)) = Ar(G'^' - G^) = -27riNv{r) , 

so that the prefactor in Eq. ([l0| ) is the product of two 
densities of states. If we now introduce the free energy 
F{uj) as 



FH = hm / VQe^pi-T[Q;co]), 

N^O l\ J 



(12) 



we see from Eqs. ( |10| ) and ( |ll| ) that R{uj) can indeed be 
written in terms of F{u)) as in Eq. (^. 

We can now apply the usual perturbative methods to 
expand F{lij) in powers of the coupling constant (inverse 
dimensionless conductance) 1/g = and then ap- 

ply the relation (||) to obtain R{lu). The great advantage 
of this method is the much smaller number of diagrams 
we need to evaluate. For example, in three-loop perturba- 
tion theory for the orthogonal case there are 5 diagrams 
for F{uj) comparedj-tp 41 for R{uj). 

In previous worked a source field for i^iE) has been in- 
troduced with a complex index structure similar to that 
necessary for calculating the conductance moments. A 
considerable simplification here is that u itself suffices 
as a source field for calculating iy{E). We note that this 
method cannot be extended to the supersymmetric sigma 
modeE where the free energy is always unity owing to 
the fact that the supersymmetry is preserved in the ef- 
fective functional. To break the supersymmetry one must 
introduce the fc-matrix in the expression for the prefac- 
tors, STr(fcA(5), or equivalently introduce an appropriate 
source term in the effective functional^. Although this is 
straightforward, all 'savings' in terms of the number of 
diagrams would be lost. Thus, for perturbative calcula- 
tions the replicated a model turns out to be more eco- 
nomical than the supersymmetric one. Finally we note 
that the free energy has previously been calculated up to 
four- loop ordeinil in order to find P{g)- Our result above 
gives a direct physical meaning to the field-theoretic free 
energy. 



To generate perturbation theory we must introduce a 
parameterization of Q that satisfies the saddle-point con- 
straints. There are several parameterizations available, 
and they will give different contributions for a given dia- 
gram, but the sum of all contributions at a given order of 
perturbation theory will always be the same. The partic- 
ular one used depends upon the application. The choice 



Q ^ A{W + ^/T+W^) 





2iV+ 



2iV 




(13) 



reproduces the results of the cpginal impurity diagram- 
matics in the diffusive regimeciJ, and is thus useful for 
direct comparison between the results of impurity dia- 
grammatics and the non-linear cr-model. Here V is an 
unconstrained N x N matrix with elements of the form 
appropriate for the given Dyson ensemble. The func- 
tional integration is then carried out over the indepen- 
dent variables V. Another parameterization,E2l 



Q = A{l + W/2){l-W/2y 



iU) 



where W has the same properties as before, is more con- 
venient for calculation because many terms then vanish, 
and the Jacobian of the transition from Q to is unity 
in the replica limit. Since the Jacobian contributes terms 
necessary to remove ultraviolet divergences, this means 
that the sum of all terms at any order in perturbation 
theory will have no ultraviolet divergence and so there is 
no need for regularization. Considerations in the ergodic 
regime are also much simpler, since Jacobian terms would 
contribute there. In what follows we will, therefore, use 
these parameterization. 

When we substitute a parameterization into the free 
energy functional J-^[Q] we obtain a sum of vertices. 



(15) 



where the J-2n\W] vertex contains 2n powers of W . The 
first term J-2 is the quadratic part, and leads to the lad- 
der propagators. The perturbation expansion up to the 
third loop order is then obtained by expanding out the 
exponential in the other J-2m 



exp {-Ti - Tq- 



and averaging with respect to exp {—T2) 
tuting the parameterization we get 



Upon substi- 



Ti^TTv j dfirTv [-2DyVVV+VV+ - iujVV+VV+] 

Tfi^TTv j d'^rTv [D{2VVyV+VV+ + yVV+VVV+VV+) - iuoVV+VV+VV+] 



(16) 
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If we then Fourier transform into momentum space we 
recover the expected diffusion propagator of Eq. (^) as 
the average of VV'^ with respect to exp (—^2), 



{V.a{q)V+{-q)) = 



1 



1 



2'KV Dq^ — iuj 



Si j Sap 



(17) 



The higher-order terms yield the effective four- and six- 
point vertices 



V4{qi) = 2711^ [-DAl - 2iuj] 



2iiuj\ 



(18) 
(19) 



where A™ is the sum of all distinct scalar products of 
pairs qi and qi+m of the n incoming momenta which arc 
separated by m vertices, 



^qi-q. 



i-\-m 



(20) 



Note that in the derivation of Eq. ( |18| ) we have sym- 
metrized over all incoming momentum variables. 

If we use the above formulas to calculate the expansion 
of free energy F{ijj) to three-loop order, the diagrams ob- 
tained are simply those we derived earlier using the tra- 
jectory language. The field theory enables us to associate 
an algebraic expression to each component in a diagram, 
and we obtain the results 



F^ = - 
F2 = 

Fja — 

Fab = 
Fzc - 
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E 



(2t:u) ^ V1V2 

91,92 

2 ^ {iujf 
{2ttvY ^ VIV2V3 

91,92,93 ^ 



1 



{2T:vf 



E 

1,92,9. 

E 



V1P2V3 

Vi - 2iuj 



(27rry)2 ^ ViV2r3 
91,92,93 



(21) 



F. 



3d 



{Pi2 - ^CJ)' 

4(27rzy)2 ^ Pir2P3Vl23 

91,92,93 



E 



Fse 



2{2Txvy 



E 



(7^12 - i<^){V\z - iuj) 
V1P2V3V123 



where we use the notation 



Vl...n = D{qi 



■ qn 



Of the above terms only Fi, F^b and -Fsd contribute in 
the unitary case. All diagrams contribute in the orthog- 
onal case, and there is an extra factor of two arising from 
the two possible relative directions of the electron lines. 

Symbolic expressions for F similar to those in Eqs. ( |2l| ) 
have been previously derived in Ref. |2^ where the param- 
eterization (|l3|) has been used. Equations (21) derived 



within the parameterization (|j) are valid for to both the 
ergodic and diffusive regimes. 

Let us first check that these formulae reproduce the 
known results for the ergodic regime to this order. The 
ergodic regime corresponds to the zero-mode approxima- 
tion in which we set all qi equal to zero. Doing this and 
differentiating twice with respect to liJ yields 



R°'\x) = 3?e 
i?""'(a;) = 5fte 



1 1 

-I- 



{ixY {ixf 2{ix)^ 
1 



2{ixf 



(22a) 
(22b) 



where x = ttw/ZV. These are the correct expansions of 
the exact resultsE Indeed, ipr K{t), which is the Fourier 
transform of R{x), one hasO in the region <t ^ 1: 

K°'*{t) = t[2 - \n{t + 1)] fii2t^t^ + (23a) 



^"'"(t) = t0{l -t) = t 



(23b) 



Here t is measured in u nits of tn = h/A. (To compare 
the results of Eqs. ( ^2h| ) and (23a), one uses the result 
that the inverse Fourier transform of t" is j"+^n!/2a;"+^.) 
Recall the famous result that there are no corrections to 
the universal RMT result in the unitary case to any order 
in perturbation theory. This happens due to the cancel- 
lation between rather than the absence of appropriate di- 
agrammatic contributions, and is therefore a useful check 
on our calculation. Note also that the complete perturba- 
tive expansion for R{uj) would give, upon Fourier trans- 
formation, the exact K{t) only for t < 1. The pertur- 
bation theory cannot give the discontinuity point, t — 1, 
because this is controlled by the non-trivial saddle pointH. 

We will now derive the leading-ordep-contributions to 
i?(w) in 2d. This was previously donellll for the orthog- 
onal case, and the form in the unitary case conjectured 
from a one-parameter scaling hypothesis. We replace the 
sums over q in Eq. by integrals and use dimensional reg- 
ularization to evaluate these integrals in d = 2 + e dimen- 
sions. (Note that all the relevant integrals may be found 
in Ref. p3). Then we carefully take the limit e ^ 0, keep- 
ing only the terms divergent in this limit. The two-loop 
orthogonal result is 



Frico) 



1 



TT^P {2ttv){2t:DY 



Taking the second derivative with respect to uj, and not- 
ing that dimensionless conductance g = Gjie^j-KK) = 
Att^i^D, gives 



i?r(^) 



2A 



(24) 



Finally we let e ^ and note that {—iujY/e becomes 
ln(— iwr). Hence 



2A 



Ke 



In (— iwr) 



1 

2 

g s 



A 
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The symplectic result would simply follow upon multi- 
plying by -1/2. 

Next we consider the three-loop unitary case. Careful 
evaluation of the integrals yields 



3A 



Ke 



In [—iujt) 



(26) 



2 rd 



1 



Finally we calculate the three- loop orthogonal result: 
A^L-^ 1 (4-3e)(-«a;)3<^/2-Hi 



Thus 



3e^ 



3e3 



so that 



5Re 



2e 



(25) 



2A (4 + 3e)(-ic^)3^/2-i 



which in the limit e ^ becomes 



7rg3 2e2 
which in the limit e — > gives 



(27) 



9A 
2^ 



3?e 



In^ (— io^r) + In {—ilut) 
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l-h21n(sAr) 



We now show that the same results can be obtained using a one-parameter scaling hypothesis in which the 
renormalized conductance, g{uj), is substituted into the lowest order R{uj) diagram. The latter has the form 
R{uj) — (A2/7r2/3)Ke 'Pi^, and substituting here D = Dq + 5D{uj) gives upon the expansion in 5D 



A2 

7r2/3 



1 



(28) 



We obtain the perturbative correction to the diffusion 
constant, 5D{ijj), from the (i{g) function which is defined 
by 



d\ng{L) 
dlnL 



To two-loop order the orthogonal and unitary beta func- 
tions are P°'^^{g) — 2/g and /3""'(5) = 2/g^ from which it 
follows that 





Effort 


Do 


50 




Sgnni 


Do 


90 



— ln(L/£) ^ 
9o 



2 {-iajy/^ 
90 e 

1 i-iujy 

5o e 



We can then substitute these results into Eq. (28), and 
use the following results for the ^-integrals. 



Doq' 



y 

^ (Dog^ - iw) 



(^09- 



^ {Doq^ 



2\2 



8ttD 

4 + 3e 
IGttD 



We find that this exactly reproduces the results of Eqs. 
(U), (^) and (HI), proving the validity of the scaling 
hypothesis to three-loop order. Note that the conjecture 
of Ref. |l^ for R™^'^{uj) does not produce the correct nu- 
merical factor in Eq. (|2^) . The reason is that the calcula- 
tion in Ref. ^ has been performed in exactly two dimen- 
sions where it was possible to take into account all weak- 
localization logarithms but not logarithmic corrections 



arising from stronger mesoscopic divergences. To pick 
up all logarithmic corrections one must work in d = 2 + e 
and take carefully the limit e ^ 0, as we did here. 

In conclusion, we have shown that the language of 
semiclassical trajectories suggests the most economical 
way of drawing diagrammatic corrections to spectral cor- 
relations in disordered electronic systems up to a high 
order of perturbation theory. We have introduced the 
weak diagonal approximation which includes the con- 
tributions of pairs of trajectories which are made from 
identical pieces joined together in different ways at some 
self- intersection (in real space) points. This gives physi- 
cal meaning to the loop expansion of the field-theoretical 
'free energy': the second derivative of this free energy 
with respect to frequency is the two-level correlation 
function. We have shown this directly by using the repli- 
cated nonlinear a model. Note that such a derivation 
does not work for the supersymmetric nonlinear a model 
as the free energy equals zero unless the supersymmetry 
is broken. Naturally, all the perturbative results may be 
reproduced within the supersymmetric model but in a 
much less economical way. Using the method described 
above, we have calculated the leading order contribu- 
tions to the two-level correlation function in 2d in the 
non-ergodic regime where the standard diagonal approx- 
imation gives vanishing results. 

In the ergodic regime we have found that the loop ex- 
pansion reproduces the deviation from linear in time be- 
havior of the spectral form factor in the orthogonal case 
known from random matrix theory. In the unitary case, 
we have demonstrated that all such corrections are mutu- 
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ally cancelled. The random matrix theory is believed to 
describe correctly level statistics of classically chaotic sys- 
tems in the universal regime. It is therefore reasonable 
to conjecture that the weak diagonal approximation in- 
troduced here should be valid for generic chaotic systems. 
Although the universal regime is well understood within 
the supersymmetric approach, the weak diagonal approx- 
imation might be extended to the non-ergodic regime in 
chaotic systems. 
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